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Abstract 

This paper is part of a program that aims to understand the connec- 
tion between the emergence of chaotic behaviour in dynamical systems in 
relation with the multi-valuedness of the solutions as functions of complex 
time t. In this work we consider a family of systems whose solutions can 
be expressed as the inversion of a single hyperelliptic integral. The asso- 
ciated Riemann surface 1Z — » C = {r} is known to be an infinitely sheeted 
covering of the complex time plane, ramified at an infinite set of points 
whose projection in the T-plane is dense. The main novelty of this paper 
is that the geometrical structure of these infinitely sheeted Riemann sur- 
faces is described in great detail, which allows to study global properties 
of the flow such as asymptotic behaviour of the solutions, periodic orbits 
and their stability or sensitive dependence on initial conditions. The re- 
sults are then compared with a numerical integration of the equations of 
motion. Following the recent approach of Calogero, the real time trajec- 
tories of the system are given by paths on 1Z that are projected to a circle 
on the complex plane r. Due to the branching of 1Z, the solutions may 
have different periods or may not be periodic at all. 

Keywords: Complex dynamics, Riemann surface, inversion of hyperelliptic 
integrals, isochronicity, sensitive dependence. 
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1 Introduction 



It has now become a classical subject to study the connections between the 
intcgrability of a dynamical system (or the absence of it) and the singularity 
structure and multi-valucdncss of its solutions. The discovery of this relation 
dates back to the work of Painleve and his collaborators [22,29] , who classified all 
second order ODEs whose solutions have movable poles as their only singularities 
(now known as the Painleve property) and are therefore single- valued functions. 
Then Kowalewskaya found a new integrable case of the heavy top by requiring 
that the solutions have the Painleve property [23,24]. A number of techniques 
collectively known as Painleve-Kowalewskaya analysis have been developed over 
the last thirty years (for a review see, for instance, [31]) to test for this property 
of the solutions, essentially by seeking for a formal solution near a singularity 
in terms of a Laurent series, introducing it in the equations and determining 
the leading orders and resonances (terms in the expansion at which arbitrary 
constants appear) . Painleve analysis has been extended to test for the presence 
algebraic branching (weak Painleve property, [30,31]) by considering a Puiseaux 
series instead of a Laurent series. These analytic techniques (which have been 
algorithmized and are now available in computer packages) constitute a useful 
tool in the investigation of integrability: in many non-linear systems where no 
solution in closed form is known, Painleve analysis provides information on the 
type of branching the general solution or special families might have. It has 
also proved to be useful to identify special values of the parameters for which 
generically chaotic systems such as Henon-Heiles or Lorenz become integrable. 

It is also natural to investigate the singularity structure of solutions of chaotic 
dynamical systems. Tabor and his collaborators initiated this study in the early 
eighties for the Lorenz system [32] and the Henon-Hcilcs Hamiltonian [18]. They 
realized that the singularities of the solutions in complex-time are important for 
the real-time evolution of the system. The complex time analytic structure was 
studied by extensions of the Painleve analysis involving the introduction of log- 
arithmic terms in the expansion — the so called ^-series — which provides a 
local representation of the solutions in the neighbourhood of a singularity in the 
chaotic regime. Their local analytic approach was complemented by numerical 
techniques developed for finding the location of the singularities in complex time 
and determining the order of branching [16]. In all the chaotic systems under 
study, they observed numerically that the singularities in complex time cluster 
on a natural boundary with self-similar structure [17]. An analytic argument to 
explain the mechanism that leads to recursive singularity clustering was given 
in [26]. Similar studies relating singularity structure, chaos and intcgrability 
have been performed by Bountis and his collaborators. Going beyond the lo- 
cal techniques described above, the emphasis was put on a global property of 
the complex solutions: whether their Ricmann surface has a finite or an infi- 
nite number of sheets. The authors of [10] proposed to classify the system as 
integrable in the first case and as non- integrable for the second case. 

These ideas actually go back to the work of Kruskal and his poly-Painlcvc 
test to detect non-integrability [25]. Combining analytical and numerical results 
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for a simple ODE, Bountis and Fokas [21] identified a chaotic behaviour with 
the property that the singularities or branch points of their solutions are dense. 
In their work, however, it is not too clear whether the criterion is that such 
points should be dense on the Riemann surface itself or that the projection 
of the branch points is dense on the complex plane. The second condition is 
obviously much weaker. 

In fact, as was indicated in [1,2], there exists a wide class of completely 
integrable Hamiltonian systems whose solutions live on an infinitely sheeted 
covering of the complex time plane with an infinite number of branch points 
that are projected onto a dense set on the complex plane and yet are isolated 
on the covering surface. 

Painleve analysis is a local technique and does not provide information on 
the global properties of the Riemann surfaces of the solutions (the number and 
location of the movable branch points, how the different sheets of the Riemann 
surface are connected together, etc.), which can be important for understanding 
of the dynamics of the system. 

Aware of this limitation, F. Calogero initiated a new line of research by 
studying the real time dynamics of a system as a path on the Riemann surface 
of its solution in complex time. In particular, Calogero and Francoise have 
shown that many isochronous systems can be written by a suitable modification 
of a large class of complex ODEs [11]. Using local analysis and numerical 
integration in two many-body systems in the plane [13,15], it was discovered 
that the periodic solutions could have a very high period, and the mechanism 
for the period increase was conjectured. However, in those systems it was not 
possible to achieve a complete description of the Riemann surface. Although it 
was proved that all singularities are finite order branch points, the presence of 
an infinite number of them could still produce aperiodic motion as the solution 
visits an infinite number of sheets. In the subsequent work [14], the Riemann 
surface of the solution was infinitely sheeted for generic values of the parameters, 
yet it could be described in full detail. Using these rather novel techniques, 
Calogero et al. were able to derive analytic expressions for the periods and prove 
sensitive dependence on initial conditions. In [14] the projection of square-root 
branch points of the infinitely sheeted Riemann surface cover densely a circle 
on the complex plane for generic values of the parameters. As a further step on 
the way to understand the connection between chaotic behaviour and analytic 
structure, it was natural to extend the study to Riemann surfaces whose branch 
points projected onto C do not cover densely a one-dimensional curve in C, but 
the whole complex plane. The candidate objects for this study are the Riemann 
surfaces associated to the inversion of hyper- elliptic integrals, first described by 
Jacobi [27]. 

This is precisely the case of the dynamical systems treated in this paper. 
The text is organized as follows: in Section [21 we introduce a family of simple, 
generally nonautonomous, dynamical systems whose solutions can be under- 
stood as lifting a circular path to a Riemann surface associated to the inversion 
of a hyperelliptic integral. Global geometrical properties of such surfaces for 
different cases are described in Section The algorithm that performs the lift 
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of the circular contour on the Riemann surface and possible behaviours of the 
corresponding solutions are described in Section Here we make a conjecture 
that all the solutions are either periodic or singular (i.e., they escape to infinity 
in finite time). 

The trajectories of the system for different initial conditions are shown in 
Section[Sl where the path on the Riemann surface is compared with a numerical 
integration of the equations of motion. Note that trajectories of complex Hamil- 
tonians have been also studied numerically by Bender et al. in [7]. Although 
the systems they study are different from the ones treated here, their approach 
is related to the one we propose and the observed complexity of the orbits is 
similar to the results of Section [SJ Finally, conclusions and directions of future 
work are outlined in Sectional 

2 Two classes of dynamical systems whose solu- 
tions travel on the Riemann surface describing 
inversion of a hyperelliptic integral 

In [20] Flashka drew attention to a simple class of integrable Hamiltonian sys- 
tems in K 4 = (x,y,p x ,p y ) with real time t, whose Hamiltonians and additional 
commuting integrals are given by the real and imaginary parts of arbitrary com- 
plex polynomials H(£,fi), £ = x + iy, /i = p x + ip y . In contrast to the case of 
the Liouville integrability, the generic real level surfaces 



are not 2-dimensional tori or products of lines and circles, but algebraic Riemann 
surfaces with some points removed. Then, instead of straight line motion on 
the tori, one has flows on such surfaces, which may have a rather complicated 
behaviour. Solutions of such systems are given by inversion of Abelian integrals 
on Te and posesses the weak Painleve property. 

Motivated by problems of mechanics and mathematical physics, we restrict 
ourselves to complex Hamiltonians 



where £ = C(t), both £ and r are complex and the prime denotes differentiation 
of ( with respect to r. 

By fixing the energy, H — E, the complex solution is formally obtained as 
the inversion of the following quadrature 



1e = {x,y,p x ,p y | H(C,fi) = E,E eC} 



(1) 




(2) 



<K 




(3) 
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where r is the initial time and £ = C( T o)- Whenever the potential V(Q 
is a polynomial or a rational function of £, the above integral is in general 
hyperelliptic, and £(t) is a multi-valued or even infitely-valued function on the 
complex plane r. 

We emphasize that a single hyperelliptic integral often appears as a quadra- 
ture in various integrable problems of mechanics and mathematical physics, 
and the last decade the problem of its inversion revokes increasing attention, 
[1,3,19,33]. 

On the other hand, by introducing real time t such that r = xt + TQ, x being 
arbitrary nonzero complex constant, the real and imaginary parts of £(r(i)) 
give a solution of the Hamiltonian system in R 4 = (x,y,p x ,p y ) with a real 
Hamiltonian proportional to Re xKeH + Imx ImH and describe flows on the 
hyperelliptic curve H(£, y) — E. 

Some properties of such flows were studied in [4] and, from the algebraic 
topology point of view, in [28] (see also references therein). 

In our paper we focus on a similar yet different class of systems that describe 
not straight-line, but circular motion on the complex plane r. 

Calogero's circular motion 

F. Calogero has shown that a large class of complex ODEs can be modified 
in such a way that the resulting system posseses a large number of periodic 
solutions [12]. Following Calogero's approach, let us construct another class of 
systems associated to the quadrature (0 by performing the following change of 
independent and dependent variables, 

r-ro = ^(e^-l), (4) 
C(r) = e— 4 z(t), (5) 

where u) > is a real parameter that plays the role of a frequency to which 
the fundamental period T — 2tt/uj is associated and a 6 1 is a parameter still 
to be determined. The effect of the transformation JSJ is to introduce the real 
time evolution as a parametrization of a circle in the complex r-plane. Thus 
the name circular motion as opposed to the straight line motion studied by 
Flaschka. If the solutions of the complex ODE were analytic or meromorphic 
functions of t, then the resulting real time evolution would be periodic, but if 
the solutions are multi-valued functions of t (possibly infinitely-valued) more 
complex behaviours are possible. 

In what follows, for the sake of simplicity, we concentrate on the case 

V(Q = C k+1 , kez+. 

The restriction of V((^) to a single power rather than an arbitrary polynomial 
is by no means essential for the analysis of the geometrical properties of the 
Riemann surface performed in Section[3] As a matter of fact, the results of this 
Section are valid for an arbitrary polynomial function V(Q . 
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Under the transformation (@J-©, the complex ODE |J2J with potential 
V(Q — C 1 transforms into 

e Ma-a)t + iw ( 2Q; - i)i + a (i - a ) w 2 z ] = + i) e i«^* ^. ( 6 ) 
If we require the transformed system to be autonomous then, necessarily, 

a=— ^— , /ceZ+. (7) 

In this case we obtain the following second order non-linear ODE for the function 

z = z{t): 

, + 1 ^_j w ,-^^^ = -( fc + l)**, (8) 

where z S C but now £ € R is an ordinary (real) time variable, and dots 
denote derivative of z(t) with respect to its argument. The complex variable 
((t) = ((r(i)) evolves in time according to the following non-autonomous ODE: 

C-i^C--(fc + l)e 2iwt C fc , (9) 

which corresponds to setting a = in J5J|. The solutions of © and JSJ are of 
course related by the simple transformation 10}, a change into a rotating frame 
of reference with constant angular velocity, and the relation between the initial 
data is: 

C(0) = z(Q), C'(0) =i(0)+ iauz(0). (10) 

If the potential V(C,) were not a simple power but a polynomial in 
Calogero's trick would lead to a system with non-autonomous t-periodic terms 
which can be interpreted as external periodic forcing. The geometrical descrip- 
tion of the Riemann surface associated to the inversion of © when V(() is an 
arbitrary polynomial of degree 5 is given in the following section. 



Three different examples. 

We will choose three different integer values of k in JHJ which lead to systems 
with quite different properties. 

1. If we set k = 2 in (JSJ) and write it out in terms of its real and imaginary 
components z = x + iy, we obtain the following system of two coupled 
second order real ODEs: 

£ + buoy — 6uj 2 x = — 3(x 2 — y 2 ) , 
y — 5ux — Qui y — —6xy. 

This system has one surprising property: it is isochronous, i.e., for any 
set of initial data {x(0), y(0), i(0), y{0)} the solution {x(t),y(t)} is a pe- 
riodic function of time with period T = — . The remarkable fact is that 
isochronicity is compatible with the non- linear character of 
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Indeed, as follows from the quadrature © for k = 2, its inversion £(r) is 
an elliptic function of r, therefore meromorphic and single- valued. Since 
r(t) is periodic in t with period T ~ 2ir/w, it follows from JSJ, Q) that 
the solution z(t) = e 2lw * C( r M) is a ls° periodic with the same period. 

2. In the second example, we set k = 4 in (JSJl which leads to the system 

x + (7w/3)y - (10lo 2 /9)x = -5 (x 4 + y 4 - 6x 2 y 2 ) , 
y - (7oj/3)x - (10w 2 /9)y = -20 xy{x 2 - y 2 ). 

In this case it is no longer true that all the solutions have the same period, 
moreover, they may be not periodic at all, since the function £(r) obtained 
by inverting the hyperelliptic integral 

is no longer single-valued. In fact (see Section^ it can be considered as a 
function on an infinitely sheeted ramified covering 1Z — * C = {r} such that 
the projections of the branch points at the complex plane form a dense 
set [27] . This example motivates the study of geometric properties of the 
Riemann surface 1Z and its covering of C. As we shall see, the evolution of 
the dynamical system (|12|) can be analyzed by following a circular contour 
r(t) lifted to the infinitely sheeted covering 1Z. 

3. For the third example we choose k = 5 in (jHJl. The equations defining the 
dynamical system are now 



x + 2 ujy - 3/Auj 2 x = -6x(x 4 + by 4 - 10x 2 y 2 ) , 
y -2ojx-3/4uj 2 y = -6y(y 4 + 5x 4 - Wx 2 y 2 ). 

and the function £(r) is obtained by inverting the integral 



(14) 



r-r = [ C - d11 . (15) 
4 y/2(E - if) 

Despite the similarity between this case and the preceding one, the be- 
haviour of the solutions is very different. As we shall see in Remark 3.1 
below, for any initial condition the solution of l|14|) is periodic with period 
either T or 2T. 

The above examples illustrate the importance of studying the detailed struc- 
ture of the Riemann surface that describes the inversion of a hyperelliptic inte- 
gral. 
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3 Inversion of a single hyperelliptic integral and 
the associated Riemann surface 



We first recall some necessary basic facts of the theory of algebraic Riemann 
surfaces. Let tk , . . . , rj^ be arbitrary distinct complex numbers. The genus 2 
surface 

r = {( v ,v)\v 2 = P 5 (ri)}, P 5 (v) = -(v-vi 1) )---(v-vi 5) ), (16) 

can be represented as 2-fold covering of the Riemann sphere C = {r/} U oo 
ramified at six branch points 

Qi = (^ 1} ,0), Q 5 = (^ 5) ,0), oo. 

Points on different C-sheets that have the same ^-coordinate have /i-coordinates 
with opposite signs. For this reason we denote the sheets as C + (upper sheet) 
and C~ (lower sheet). The passage from one C-sheet to the other is realized 
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Figure 1: The two sheets of T5 with the branch points and cuts (solid lines). 

through any of the three cuts that join the pairs (Q2, Qs), {Q3, Q4) and (Qi, 00) 
on C + and C _ and that do not intersect each other. Let us fix a canonical basis 
of cycles 01, a 2 , bi,b 2 on T such that 

<Zj o cij — hi o bj =0, a,o bj = 6 i3 , i,j = 1, 2, 

where 71 o 72 denotes the intersection index of the cycles 71,72- Following the 
common convention, a-cycles lie entirely on the upper sheet C + and embrace 
one of the cuts while 6-cycles run through both sheets. We assume that a\ 
embraces (Qi,oo) and b\ joins 00 with Q3, so that the projections of cycles 
ai,&i embrace 00. For the pentagonal curve T$ = {{rj.^)\ii 2 = r/ 5 — E}, the 
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77-coordinates of the finite branch points are 

^ ) =exp(^fcH + iarg £ ;)^ fc = l,...,5. 

For real E > 0, the branch points and cuts are depicted in Figure 1, where 
are shown as planes and Qj are denoted just by their indices j. Figure 
2 also shows the canonical cycles, where segments of the cycles 61,62 on the 
upper and lower sheets are given by solid and dashed lines respectively. Since 
the branch points (Qi, 0), i = 1, . . . , 5 are just the five complex roots E 1 / 5 , for 
generic (complex) values of E their position is simply obtained by a rotation 
and scaling of the pentagon of Figure 2. 
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Figure 2: Canonical cycles a, and bi, i — 1,2 for the pentagonal curve T$ = 
{(77, fi) \ /i 2 = i] 5 — E}. Segments of cycles on the upper and lower sheet are 
depicted as solid and dashed lines respectively. 



We introduce next the following basis of holomorphic differentials 




on r and consider in detail the problem of inverting the hyperelliptic integral 

r=/ P ^-, P = ( w )er. (17) 

Joe y/2fi 

The integral has four independent complex periods defined by: 

V 2 =I^-, V S =I^-, %=/-^-. (18) 

J ai V2/Lt Ja 2 V2/x Jh x V2/J Jb 2 V2/i 

It follows that the inverse complex functions tj(t), h(t) should have four gen- 
erally independent periods on C. This is not possible, a fact that was already 
noted by Jacobi (see [27] for an historical account). If the periods V\, . . . , V4 
are not commensurable, then the functions rj(r), fi(r) cannot be singled valued 
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on C = {t}, but they can be regarded as single- valued on an infinitely sheeted 
covering 72. i — * C. This Ricmann surface has an infinite number of ramification 
points whose projections onto the r-plane form a dense set. 

Remark 1 According to the Weierstrass-Poincare theory of reductions of Abelian 
functions, (see e.g. the surveys [5,6]), the periods ofdrj/fi may become com- 
mensurable in special cases when the genus 2 curve T covers an elliptic curve £ 
and the differential drj/fi is a pull-back of a holomorphic differential on £ . Then 
1Z i— > C is just a finite order covering. For genus 2 curves admitting a finite or- 
der automorphism, a complete list of their Riemann matrices that characterize 
such coverings was given by Bolza in [9] (see also [8]). 

We note that the pentagonal curve T§ = {ji 2 = E — 77 s } does not belong 
to this list, hence its periods are incommensurable ( see also Proposition ^ and 
relation kilty) below). On the contrary, the sextic curve Tg = {/i 2 = E — r] e } does 
appear in Bolza 's list: it is a 2-fold ramified cover of two elliptic curves, 

£\ < r 6 > £ 2 , 

£ 1 = {w 2 = E- z 3 }, £ 2 = {wl = z 2 (E - z 3 )}, 

(Z 1 ,W 1 ) = 7Tl(?7,/i) = (77 2 ,/i), (z 2 ,W 2 ) = 7T 2 (77,^) = (l] 2 ,^^). 

Then drj/fi is the pullback of the holomorphic differential dz 2 /w 2 on £ 2 . As a 
result, the corresponding covering TZ 1— > C is only 2-fold, which implies that, 
similarly to the case of the system all the solutions of equations are 

periodic, either with period T = 2tt/lu or 2T. 

Remark 2 Note that for an arbitrary curve T given by $16p . the differential 
dn/n has a double zero at 00 and no zeros elsewhere. Namely, choosing a local 
coordinate S(P) in a neighborhood of 00 G T such that 5 (00) — and r\ = 1/d 2 , 
= 0(1/S 5 ), one has drj/^i = 0(S 2 ) dd and, therefore, t = 0(S 3 ), 6 = {\/t). 
This implies that the functions 77(t),/^(t) have only critical points r» with the 
behavior 

n(r) = (r- t»)- 2 / 3 ( Co + Cl (r - r») 2 / 3 + c 2 {r r») 4 / 3 + •••), 
MM = (r - T,)- 5 / 3 (d + di(r - r,) 2 / 3 + d 2 (r - r,) 4 / 3 + ■ • • ), 

where Cj and dj are constants. Thus, the singularities t* are 3 rd order branch 
points of the covering 72. f — *■ C and there is no other type of branching. 

Despite the complexity of the Riemann surface 1Z, it is possible to give an 
explicit description of its geometric and topological structure, which enables 
us to describe the behavior of the solutions of system (|12H . Such an explicit 
description of 1Z is new to our best knowledge. 

Description of a sheet of 1Z 

We start by making some preliminary observations. Namely, on the C + -sheet of 
r we make two extra cuts and B 2 that join 00 with Q\ and Q3, respectively 
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(see Figure|3J). We assume that B^~ and B^ are embraced by projections of the 
cycles a\ and b\ . It is always possible to choose the cuts in such a way that path 
integral l|17fl along each of them yields a straight line segment on the complex 
r-plane 1 . We also introduce the two cuts B^ and B^ on the C~-sheet with the 
same properties (in Figure |3| such cuts are depicted as the bold dotted lines.) 




Figure 3: Canonical cycles a* and bi and the new cuts B\ and B2 (dashed-dotted 
lines) needed to define the B-rule. 



Definition 3 (B-rule) Consider the hyperelliptic integral |J7|) and allow the 
point P to range over the whole curve T in such a way that the integration path 
from 00 to P does not cross any of the cuts Bf 1 and B^ . We shall refer to this 
condition as the Border rule, or for simplicity, the B-rule. 

Note that since the cycles ai and b\ cross the cuts B2 and B\ respectively, 
the corresponding periods V\ and V3 are not allowed by the B-rule, and we shall 
refer to them as the forbidden periods. Consequently, we shall refer to the other 
periods V2 and V4 as the allowed periods. 

Proposition 4 Under the B-rule, the values of the integral t{P) range over 
the whole complex plane C except for an infinite number of non-intersecting 
identical domains Wij, i,j G Z which shall be called windows. Each window is 
a parallelogram spanned by the forbidden periods V± and V3, while the window 
lattice Wij is generated by the allowed periods V% and V4 . The corners of the 
windows correspond to the image of the point at 00 € T. At each pair of points 
Ti , T2 lying on the opposite edges of Wij such that T2 = t\ + V± or T2 = t\ + V3 
the coordinates 77, [i take the same values. 

1 Note that generally the integration along a straight line segment on the rj-plane produces 
a curvilinear segment on the T-plane. The only exception is the integration along a segment 
on the real axis that does not contain the branch points of T. 
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In other words, the window Wij+i (Wi+ij ) is obtained from Wij shifting 
by period V% (V4). When PeT makes a linear combination of cycles 02, 62, the 
image r(P) shifts by the corresponding combination of V2, V4. More specifically, 
the two edges of the windows parallel to V\ (V3) are the image of the two opposite 
sides of the cut Bf {Bf). 

It follows that under the border rule, the two inverse functions T)(t) and h(t) 
are quasi- elliptic, in the sense that they are doubly periodic with the allowed 
periods V2 and V4, but they are not defined in the whole complex r-plane. The 
domain of definition T> = C \ {U Wij} can be regarded as one of the infinite 
set of sheets of the Riemann surface 1Z. Inside each parallelogram of periods V2 
and V4 the functions have 4 poles lying at the corners of a window Wij . This is 
the equivalent of a fundamental domain for the "quasi-elliptic" functions rj(r) 
and /i(r)under the B-rule. For a general genus 2 curve T with two real periods, 
the domain T> is depicted in Figure 0] as the grey shaded area. 




r(Q.) V 2 y, 




Figure 4: Domain V (shaded area) and windows Wij. We can consider this to 
be one of the infinite number of sheets of the Riemann surface 1Z 

Proof of Proposition |4j Let us first study the integral t(P) when P e T 
moves along the cuts B^ and Bf without crossing them. When P moves on 
C + from 00 to Qi along the upper edge of B^ and returns to 00 along the lower 
edge of 2?i~ , the point t(P) performs a straight line motion along the forbidden 
period vector V±. If P continues moving on C + from 00 to Q\ along the lower 
edge of B2 and returns to 00 along the upper edge of B^ on C~, the integral 
r(P) traces the second forbidden period vector V3. If P traces the same two 
motions now on the opposite C-sheets of T, the point r(P) moves subsequently 
along the vectors — V\ and — V3. Corresponding to this sequence of motions of 
t along the cuts the integral r(P) traces a parallelogram W spanned by the 
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forbidden periods on the r-plane. The coordinates rj, fi of P coincide along the 
points on the upper edge of B± (B± ) identified with the the lower edge of B-f 
(Bf), and they also coincide along the opposite points on the upper and lower 
edges of B^ (B^). Under the B-rule the point t(P) cannot reach inside of a 
window W. Outside of W, r(P) behaves as an elliptic integral of the first kind 
with periods V2 and V4, therefore t(P) cannot reach inside of an infinite number 
of windows Wij obtained from W by translations by the allowed periods. Since 
the map r(P) is analytical, in C = {t} each Wij has a small neighborhood 
covered by the image of P, hence the parallelograms cannot intersect. Finally, 
the proof that t(P) reaches any point of T> = C \ {U Wij} goes along the same 
lines as that for any elliptic integral of the first kind. □ 

3.1 Gluing the P-sheets together. 

Let us first introduce some notation for the different edges of a window Wij 
on V: we shall label them as N ("North"), W ("West"), S ("South"), and E 
("East"). By convention, the N- and S-edges are spanned by the period Vi, 
while the W- and i?-edges lie along V3. Next, we assume that the B-rule can 
be violated and we allow the integration path from 00 to P to cross one of the 
cuts B^,Bf- only once. If the path crosses -B-f, (B^), then the point leaves the 
domain V through S-edge (respectively, TV-edge) of a window Wij . 

In view of Proposition 0] this motion can be viewed as the passage to a 
different sheet V of the surface TZ, which is a translate of an identical copy of 

V and has the same window lattice structure. On crossing the S-edge of Wij 
on T>, the image t(P) immediately arrives at the TV-edge of a window W(j of 

V and then continues travelling on this new sheet. That is, the iV-edge of a 
window W(j on V is glued to the S'-edge of Wij on T> in such a way that the 
corners are connected to corners and at the identified points the values of 77, /j 
coincide. Next, on crossing the iV-edge of the window Wij on T>, t(P) enters 
another sheet T>" from the S'-edge of a window W"j on T>" . Similarly, when 
the integration path on T crosses -B^", (B2), then t(P) leaves V through the 
E- (respectively, W-) edge of the window Wij and arrives at other sheets of TZ 
from the W- (respectively, E-) edges of windows on these new sheets. 

Since there is an infinite number of windows on V, one might think that 
each sheet is connected to an infinite set of different 2?-sheets of TZ. However, 
the actual structure of TZ is simpler and completely described by the following 
short proposition. 

Proposition 5 The Riemann surface TZ is a union of an infinite number of 
identical T>-sheets glued to each other along the corresponding opposite edges of 
their windows. Any arbitrarily chosen "home " sheet T>h is connected to precisely 
four other sheets T>hn,'E > he,'DhSi'Dhw md the passage from T>h to them is 
realized by crossing respectively the S-, W-, N-, and, E-edge of any window 
Wij on T>h- 

Proof of Proposition. We only need to prove that passage through S— (or 
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W- or N- or E-) edges of different windows of T>h leads to the same P-sheet of 

n. 

Let pit) be a path on V that starts at a point Pq , crosses the border Bf- U 
once and ends at a point P\ . Next, let p^ (t) be the path p(t) followed by a linear 
combination of cycles a 2 , &2 and p^ (t) be the same linear combination of cycles 
(with starting point P ) followed by the path pit). Then the integrals r{p^ l \t)), 
T(pW (t)) give paths 71, 72 on 1Z that start at the same point on the sheet T>h and 
leave it through the same edges of different windows on T>h ■ By construction of 
p^\t) andp( 2 )(t), the endpoints of cycles 71 and 72 have the same r-coordinate, 
as well as the same 77, /^-coordinates (corresponding to Pi e T). Hence, the 
endpoints on 1Z coincide, and by continuity the corresponding 2?-sheets can be 
identified. □ 




Figure 5: Window lattice on T>h (bold) and on T>hn (dotted). The path 
leaves T>h through the S'-side of a window (bold arc) and enters T>hn , then it 
continues on T>hn (dotted arc) and leaves the sheet through the TV-edge of a 
window, coming back to T>h 

As a corollary of the above proposition, we note that each of the four T>- 
sheets adjacent to T>h is in turn connected in the same way to other four sheets, 
one of them being the initial sheet T>h- In particular, when t(P) leaves T>h 
through the S-edge of a window Wy, enters the sheet T>h,n and then leaves it 
through the TV-edge of another window in T>h,n, then it returns to T>h (see 
the illustration in Figure |SJ|. It is then natural to label the P-sheets connected 
to T>hn as Vhnn, 'Dhne, 'Dhnw, and Vhns = T^h- 

Finally, we also state the important 

Lemma 6 The passages between T> -sheets in the "horizontal" (E, W )-direction 
and the "vertical" (N , S )-direction commute. That is, the following sheets 
coincide: 

T>HNE = T>HEN, T>HNW = T>HWN, T>HSE = T>HES, T>HSW = ^HWS- 
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Proof of Lemma. 

Indeed, the windows on T>hne and on T>hen are obtained from those on T>h 
by shift by the same vector V\ — V3, hence T>hne and T>hen can be identified. 
Similar arguments apply to identify the other pairs of sheets. □ 

As a result, the Riemann surface TZ can be described as a Z 2 -countable set 
of ©-sheets, that is 

ft = {u£>y I i,j e z} 

where i,j enumerate the passages in the "horizontal" and "vertical" directions 
respectively. Hence, travelling on TZ can be encoded as a path (sequence) on a 
Z 2 -lattice A. 

3.2 Branching of the projection 1Z C = {t}. 

Note that each corner of a window on a 2?-sheet is connected to several different 
sheets, hence the corners are branch points of the projection II : TZ 1— > C = {r}. 
In particular, consider the corner between the E- and S'-edges of W%j on T>h- 
Let 7* be a path on TZ that starts on T>h and whose projection to {t} is a 
sufficiently small circle around t* with a counter-clockwise orientation. Then 
7* leaves T>h through the S'-side of a window and visits the sequence of sheets 

T>H — > T>HN — * 'DhNe{= DhEN) — ► T^HE — > 

Thus r* connects four 2?-sheets and four edges of windows on them, and the path 
7* closes after making 3 complete rotations about r*. This gives a geometric 
interpretation of the cubic root branching of functions t](t) and /i(r) in (|19|) . 
The other window corners have the same properties. Note that there is no 
contradiction between the fact that each branchpoint t* is of the third order 
but connects four 2?-sheets together. The reason is that the 25-sheets are not 
copies of the whole complex plane C, and due to the presence of the windows it 
is not possible to turn by 2tt around T* on a D-sheet. 

Due to the incommensurability of the size of the windows (the modulus of 
the forbidden periods \V\\ and | V3) with the shift vectors that generate the 
windows lattice (the allowed periods V2 and V4), the projection of the corners 
of all the windows in all the sheets of TZ (i.e., the branch points of II) is a dense 
set in C. However, as we have seen above, on each of the infinite P-sheets of TZ 
the branch points are isolated and they form a regular lattice. 

3.3 Duality. 

The Z 2 -lattice A introduced above has a much richer structure. More precisely, 
one can extend it to a cell structure A on R 2 : the adjacent vertices of the lattice 
are regarded as vertices of square cells A fc; on M 2 . 

There is a remarkable duality between the objects of TZ and A: a 2-dimensional 
P-sheet in TZ corresponds to a vertex of A, the edges of the windows connecting 
two different P-sheets of TZ correspond to the segments joining the two corre- 
sponding vertices. Finally, a branchpoint r* of TZ which connects four P-sheets 
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corresponds to the square cell A^z in A defined by the four vertices which cor- 
respond to the sheets. 

3.4 The surface 1Z as the universal covering of the theta- 
divisor. 

Now, instead of the single hyperelliptic integral (|17|l . consider the vector of 
holomorphic integrals 

i \T f P ( dr l vdr]\ T 

(ui,u 2 ) = / — , — , Per, 

which has 4 independent vectors of periods in C 2 (ui,U2) with respect to the 
cycles oi, . . . , &2- As it is known from the theory of Abelian varieties, the map 
A : P i—* {u\,U2) T describes a well defined smooth embedding of T into the 
Jacobian variety Jac(r) = C 2 /A, where A is the lattice generated by the above 
period vectors. The image of the embedding is known to be a one-dimensional 
analytic subvariety 0' C Jac(T), a translate of the theta-divisor defined as the 
zero locus of the Riemann theta-function 9{u\, u%) of the curve Y. 

Let C C 2 be the universal covering of 0. Note that it can also be re- 
garded as an infinitely sheeted ramified covering of C = {u{\ and that the first 
coordinate u\ jointly with the point P e V define a point of uniquely. It 
follows that by identifying u± with the complex coordinate t in (|17f) . one can 
also identify the universal covering with our Riemann surface 1Z: two points 
on both complex surfaces are identified if and only if they are projected to the 
same u\ — t and correspond to the same point on Y . 
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3.5 Periods of the pentagonal curve. 

The existence of a group of automorphisms of the curve r 5 = {p 2 = E — rf} 
implies some special relations between its periods, which we did not find in the 
literature and indicate here. 

Proposition 7 For the chosen basis of cycles on T5, the local dependence of 
the periods on the complex energy is as follows 

V l {E)=E- z/w V? 1 i = l,..., 4, (20) 

where V® are the periods for E — 1 given by 

V? = 2 [ - dv = 1.819915 i 
Joe V^rf 

o _ 3 + V5 
V2 - —Vn (21) 

_ 3 + V5 Q 
K* - 2 V3 , 

Note that when £J varies along a loop embracing the origin and returns to 
its original value, the curve T$ undergoes the automorphism p that rotates its 
finite branch points by angle 2ir/5, whereas the periods undergo a monodromy: 
they are all multiplied by exp(— Sni/5) and, on the other hand, become linear 
combinations of V®, ■ ■ ■ , V4 3 with integer coefficients. Proposition (|21|l implies 
that for any complex energy E the following relations hold 

t = t = -^T^ Wi\ = Ws\, W 2 \ = \V 4 \. (22) 

The two forbidden (allowed) periods have the same length, moreover the for- 
bidden period V\ (V3) on C is proportional to the allowed period V% (V4). Since 
the proportionality constant is an irrational number, the period vectors are not 
commensurable, as was also claimed in Remark 2.1. 

It is worth noting that the pentagonal curve T§ is rather special, and the 
form of the windows on T> for a generic genus 2 curve {p 2 — P^irj) — E} depends 
on £ in a nontrivial way: Vj form a fundamental basis of solutions of a fourth 
order Picard-Fuchs ODE with the independent variable E. 

A sketch of the proof of Proposition [7| 

Under the subsitution 77 = i? 1 / 5 A, the holomorphic differential drj/fi takes 
the form E~ 3 / 10 dX/y/1 — A 5 , while on the A-plane the branch points and the 
integration cycles are the same as on ?7-plane for E — 1. This implies the 
relations J20J- The calculation of the periods V®, ■ ■ ■ , V® follow the ideas of [9]. 
More specifically, the monodromy transformation described above implies that 
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the periods V/ 1 must satisfy a system of homogeneous linear equations with 
integer coefficients 



Expressing the deformed cycles p(a\), . . . , p(&2) m terms of the original cycles 
on r$ for E = 1, it is possible to calculate the coefficients an,. Finally, solving 
this system and numerical integration of V® yields expressions l|21() . □ 

4 Circular Travelling on the Surface 1Z 

As mentioned in Section 1, one can consider a complex solution of the system 
(1.12) with t € K as a smooth path j(t) on the infinitely sheeted surface 1Z 
corresponding to the curve T such that the projection H(j(t)) is the circle C on 
the r-plane parameterized by (0J and thus having radius 1/u. 

4.1 A geometrical algorithm 

Comparing the quadrature l|13|l and the integral (|17f) . we find that for initial 
complex conditions Co, Co, the T-coordinate of the starting point 70 = 7(0) is 



where the energy E is defined in and the sign of the root is chosen in 
accordance with the sign of C'(0). 

Geometrically, 70 can be chosen on any 2?-sheet by adding To to a corner of 
any window 2 . By convention, the chosen initial sheet is labelled as 2?oo- 

Let Co be the preimage of the circle C on T>oo- Naturally, if Co does not 
cross any window on 2?oo, then the path II -1 C on 1Z remains on 2?oo and the 
solution C(i) of © is periodic with the fundamental period T = 2-k/uj. In other 
cases the path j(t) gives rise to a finite or infinite sequence on the Z 2 -lattice 
{UPy I i,j £ Z} according to the following algorithm: 

1) Given position of the circle Cq on 2?oo with the initial point 70 € Co, one de- 

termines the points of its intersection with the edges of Wy . Among them 
one chooses the point 71, which is closest to 70 in the counterclockwise 
direction on Co- 

2) Given 71 and the edge containing it, one finds the point j[ on the opposite 

edge of the same window such that 

7i = 7i ± Vi or 7( = 71 ± V 3 . 

Then j(t) passes to the new sheet V = 2?o,±i (respectively, V = 2?±i,o) 
with 7J as the starting point on T>' . 

2 If for a chosen corner to appears inside a window, then to should be added to another 
corner. 



exp(-37ri/5)V^° = anVi H + anVi, 



Oik G z - 




(23) 
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3) On the new sheet V one constructs the circle C\ by translating C by ±V\ 

(respectively ±14) so that j[ € C\. One then determines the intersections 
of C\ with edges of windows on T>' and chooses the intersection point 72, 
which is closest to 7J in the clockwise direction on C\ . 

Then the procedure repeats as in step 2) with 2?oo & n d 71 replaced by V 
and 72. 

4) If after n iterations one returns to 2?oo, then the corresponding circle C„_i 

necessarily coincides with Cq. If the segment of Co between j n and 7 n +i 
contains 70, then the procedure stops and j(t) is closed. Otherwise, the 
algorithm continues on step 1) with 70 replaced by j n . 

This geometric algorithm for travelling on the fundamental domain includes 
the essential properties of the flow defined by the system l|12J) and allows to 
predict the period of the orbit without numerical integration. 

4.2 Motion on a fundamental domain T>. 

There is an equivalent description of the procedure. Namely, let T> be a funda- 
mental parallelogram of allowed periods V2 and V4 on T> that contains a window 
W* completely inside. Then T> does not have any other windows inside and 
there is a natural bijection between the fundamental domain T> = T> \ W* C T> 
and the curve V. A smooth path j(t) = II -1 C on 1Z induces a smooth path P(t) 
on r and a generally discontinuous path ^(t) on T> that consists of segments 
of translates of Cq: when ^(t) crosses an exterior edge of 2? at a point 7*, it 
reappears on the point 7** on the opposite exterior edge of the domain, since 
both points correspond to the same point on T. Similarly, on crossing an edge 
of the interior window of T>, the trajectory continues from the opposite edge of 
the window as described in step 2) of the algorithm. 

Example 

For the pentagonal curve with the complex energy E = 1110 + i 1452, a 
fragment of the initial sheet 2?oo and the corresponding fundamental domain V 
are presented in FigureEK), and^p) respectively. Here we also show an example 
of the circle Co on 2?oO; which generates a periodic path ~/(t) in T> with period 
IT . The corresponding periodic sequence on the Z 2 -lattice is 

(o^) - (-1, 0) (0, 0) -> FiTo) -> (-1, 1) -> (071) (0, 2) -> 

(13)^(1,1)^(1,2)^(1,1)^(2^)^(1,1)^(2,1)^ 
(2,0) _(o^i)-> (0,0)- (0,-1) -(0^)- 

where the bars mark the 2?-sheets visited when the projection ITy(£) performs 
a complete turn. Applying the above algorithm to different initial conditions, 
we obtained a great variety of periodic trajectories. In all cases we observed a 
perfect coincidence of the motion in T> with the results obtained by the numeric 
integration of the equations © and \V2\. 
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Figure 7: (a) Circle Co on the initial sheet 2?oo- (b) Path j(t) on the funda- 
mental domain D. (c) Sequence of sheets of 1Z visited by "fit), (d) Orbit of 
obtained by numerical integration of @ with period 7T. 

4.3 Aperiodic and isochronous trajectories 

The algorithm described above does not cover the cases when, at a certain step, 
7(i) arrives precisely at the corner of a window. Since the corners are branch 
points of 7Z connecting 4 different ©-sheets, it is not possible to choose uniquely 
the next sheet, and at this step the procedure should terminate. Note that in 
this case the corresponding complex solution z(t) escapes to infinity in finite 
time. 

We also observe that if an aperiodic path j(t) exists, then it must pass 
arbitrarily close to some branch point on 1Z for sufficiently large time. The 
following obvious properties hold. 

Proposition 8 If the path 7(i) = If -1 C is not periodic, then it visits an infinite 
set of T> -sheets, and for any S > there exist a time t* and a branch point t* € 1Z 
such that \j(t*) -t*\ < 5. 

In other words, if an aperiodic orbit exists, then it has sensitive dependence 
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on initial conditions: if 70 yields a non-periodic path j(t) on 1Z and its arbitrarily 
small variation 70 results in a path j(t), then there always exist a finite t* such 
that for t < t* both paths are very close to each other, but for t > t* they pass 
by a branch point r* on different sides. As a consequence 7(4) and j(t) enter 
different P-sheets and their future behavior is completely different. 

Proof of Proposition 151 

Suppose that there exists a non-periodic path j(t) visiting a finite number 
of 2?-sheets, then some of these sheets are passed an infinite number of times. 
On each of them the corresponding translate of C is the same and consists of a 
finite number of arcs. Hence at least one of the arcs is passed more than once 
and therefore j(t) must be periodic. For the second part, since j(t) visits an 
infinite number of 2?-sheets, its image in the fundamental domain T> consists 
of an infinite set of segments of translates of C by non-commensurable periods 
V\, . . . , V4. Hence the arcs cover a dense set on D, and in particular they pass 
arbitrarily close to any corner of the window in T>. □ 

On the contrary, it is easy to prove that the periodic solutions of are 
isochronous, i.e. , for any initial data £(0),£'(0) that give rise to a periodic 
orbit, there exists an open set of initial data that contains C(0), C'(0) such that 
all trajectories starting on it are also periodic, with the same period. This result 
follows simply from the fact that branch points are isolated on 1Z and each sheet 
is connected to a finite number of sheets. 

Although one cannot exclude that non-periodic trajectories of the © and 
<|12|) may exist, we failed to find any concrete scenario of an aperiodic motion 
on Z 2 except the case of critical values of E, for which the curve T becomes 
singular. We formulate thus the following conjecture. 

Conjecture 9 For any non-critical energy E the system 111 2(1 possesses only 
periodic trajectories or singular ones ( that escape to infinity in finite time ) and 
almost all of them are periodic. However, for any k £ N, initial data can be 
found such that the trajectory has period larger than kT , where T — jio is the 
fundamental period. 

5 Numerical Results 

In this Section we want to test the validity of our results by comparing them 
with a numerical integration of the equations of motion. 

5.1 Results for k = 4 

We study trajectories of dynamical system 112f) corresponding to different ini- 
tial data. For each set of initial data {x(0), y(0), x(0), jf(0)} we compute the 
"energy" E using (|10|l and 0J. The value of E determines uniquely the four 
period vectors, which are defined generically through the contour integrals (JTSJ. 
For the special case of the pentagonal curve the four periods vectors are simply 
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given by the formula (|21[1 . As explained in Section [31 the allowed periods V2 
and V4 define a quasi- elliptic period lattice on the plane while the forbidden 
periods V\ and V3 define the windows. In the special case of the pentagonal 
curve l|22|) applies, so that each sheet of the Riemann surface TZ can be viewed 
as a fundamental cell with a window such as those depicted in Figure |5~T1 

In order to apply the algorithmic procedure described in Section 0] to lift 
the circular path to the Riemann surface TZ it is only necessary to calculate the 
position of the circle relative to the windows (the radius of the circle is always 
fixed for every initial data, and equal to 1/w), that is to say, calculating To in 
|@J by performing the integral (|23[l . 

We have considered four different sets of initial data. In all of them the initial 
velocities i(0) = y(0) — 0, while the initial positions {x(0),y(0)} are displayed 
in Table ^ For simplicity we have also set u> = 2tt so that the fundamental 
period T = 1. 



Table 1: Different sets of initial data for k = 4 



Case 


Initial data (xo,yo) 


E 


TO 


Period 


1A 


(1,2) 


67.318- 73.092i 


0.179+ 1.242i 


1 


IB 


(2,1.5) 


-112.73- 60.044i 


-0.439+ 1.09H 


3 


1C 


(1.8,2.1) 


-52.998- 215.3H 


-0.141 + 0.988i 


8 


ID 


(5.65,2.1) 


-1888.3 + 7576.3i 


0.171 + 0.140i 
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For each of the values above, the evolution circle on the window lattice of 
the home sheet 2?oo has been displayed in Figure [SU - Note that the size of the 
windows relative to the circle decreases as \E\ increases, see (|20"1) . 

The algorithm described in Section 0] to follow the path on the Riemann 
surface TZ has been implemented in Mathematica. For each of the cases 1A-1D 
the paths on TZ are all periodic, and the periods predicted by the algorithm are 
indicated in the last column of Tabled hi Figure lS~Tl we show the discontinuous 
path 7(i) on the fundamental domain T> for each of these cases. 

The path in Figure 15.11 is discontinuous because many of the arcs are con- 
tained in other sheets T>ij of TZ. Let us consider a sequence of elements of Z 2 
corresponding to the index of the sheet T>ij visited by the path after every com- 
plete turn. For each of the cases 1A-1D of Table this sequence is depicted in 
the first column of Figure 15.11 Note that the sequence in Z 2 is periodic if and 
only if the orbit of the systems of ODEs (|12|) is closed. 

Finally, we have performed a numerical integration of the equations of motion 
<|12l) for each of the initial data in Table 2] using an embedded Runge-Kutta 
method of order 5(4) with automatic step size control, developed by Prince and 
Dormand in 1981. The results of the integration show periodic orbits whose 
periods are in perfect agreement with those predicted by the algorithm of Section 

On the right column of Figure l5.1l we have displayed the orbits of the function 
((t) whose evolution is described by © and is trivially related to z{t) = x{t) + 
iy(t) by ©. 
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Figure 8: Evolution circle projected on the home sheet 2?oo f° r different initial 
data gathered in Table . 



5.2 Results for k = 5 

We study now trajectories of the dynamical system 114|) corresponding to differ- 
ent initial data. As mentioned in Remark 3.1, this dynamical system describes 
circular travelling on a Ricmann surface 1Z with only two sheets. Therefore al- 
most all solutions of Ijl4(l are periodic with period either T or 2T. This predic- 
tion is confirmed by the numerical integration of the system and some examples 
of orbits for different initial data have been gathered in Table [21 I n ah four 
cases 2A-2D the initial velocities i(0) = y(Q) = 0, while the initial positions 
{x(0), y(0)} are displayed in Tabled For simplicity we have also set ui = 2ir so 
that the fundamental period T = 1. 



Table 2: Different sets of initial data for k = 5 



Case 


Initial data (xq, yo) 


\E\ 


Period 


2A 


(1.8,1.9) 


3.350- 


10 2 


1 


2B 


(2.8,2.5) 


2.865- 


10 3 


1 


2C 


(4,4.5) 


4.781 • 


10 4 


1 


2D 


(5,8) 


7.052 • 


10 5 


1 
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Figure 9: Path 7(i) on the fundamental domain for different initial data gathered 
in Tabled 



Although all the orbits 2A-2D have the same period, a look at Figures ITT1 
and If 21 reveals that the complexity of the orbit increases with increasing energy. 

6 Conclusions 

In this paper we introduced a class of dynamical systems on ]R 4 whose solutions 
can be interpreted as lifting a circular path on the complex plane to a generally 
infititely sheeted Riemann surface TZ describing the inversion of a hyperelliptic 
integral. For the integral of the holomorphic differential di]/fi on a generic 
genus 2 surface T, we provided an explicit description of the global structure 
of TZ and showed that the surface can be identified with the universal covering 
of the theta divisor of Jac(r). (It should be stressed that in case of another 
holomorphic differential, in particular XdX/fi, the structure of TZ and type of 
branching TZ i— > C are different.) We also noted that for some regular curves 
with automorphisms, like the sextic curve Tq, the four periods of the differential 
are commensurable, and the corresponding surface TZ becomes finitely sheeted. 
This implies that almost all the trajectories of the corresponding systems are 
periodic with period T or 2T, although the trajectories increase in complexity 
as \E\ increases. 

Using the global description of TZ, we provided a geometric algorithm for 
the circular travelling on TZ, that replaces the numeric integration of the corre- 
sponding ODEs and associates to the path 7(4) on TZ a sequence on a Z 2 -lattice. 
We conjecture that for a regular value of the energy, all solutions are periodic, 
although the periods can be an arbitrarily high integer multiple of the funda- 
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mental period. 

For didactic reasons we have restricted to genus 2 curves T, although it is also 
possible to describe the global geometric structure of Riemann surfaces 1Z for 
holomorphic integrals on a generic hyperelliptic curve of genus g > 2. Namely, 
dividing its 2g periods into two allowed and 2g— 2 forbidden ones, one constructs 
an analog of a 2?-sheet by removing from C the union of identical (4g — 4)-gonal 
windows {Wij | i,j G Z 2 } which have pairs of parallel edges formed by the 
forbidden periods. Different windows are obtained from each other by shifts by 
the allowed periods. Then one can show that the surface 1Z is itself a union of 
1? 9 ~ 2 identical copies of T> glued to each other along the opposite edges of the 
windows. A trajectory j(t) on TZ can be encoded as a path on a Z 2s ~ 2 -lattice. 
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Figure 10: Left column: sequences of greets of 1Z visited every complete turn. 
Right column: orbits of ((t) obtained by numerical integration of the system 
|J5J with k = 4 for different initial data gathered in Tabled 




Figure 11: Orbits of £(t) obtained by numerical integration of the system 
with k = 5 for different initial data gathered in Tabled 
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2A 




Figure 12: Time evolution of Im£(i) obtained by numerical integration of the 
system © with k = 5 for different initial data gathered in TableEl All solutions 
are periodic with the same period T = 1 but their complexity increases with 
increasing energy \E\. 
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